Investigating the pre-main sequence magnetic 
chemically peculiar system HD 72106 



Folsom, C. P. 1,2 , Wade, G. A. 1 , Hanes, D. A. 2 , Catala, C. 3 , Alecian, E. 3 , Bagnulo, S. 4 , 
Boehm, T. 5 , Bouret, J.-C. 6 , Donati, J.-F. 5 , Landstreet, J. D. 7 

1 Dept. of Physics, Royal Military College of Canada, PO Box 17000, Stn Forces, Kingston, Ontario, Canada, 
^6"; K7K 7B4 

2 Dept. of Physics, Engineering Physics & Astronomy, Queen's University, Kingston, Ontario, Canada, K7L 
3N6 



O 

(N 
O 

Q 



O 



X 



3 



Observatoire de Paris, LESIA, 5 place Jules Janssen, 92190 Meudon CEDEX 



4 European Southern Observatory, Casilla 19001, Santiago 19, Chile 



5 Observatoire de Midi-Pyrenes, 14 Av. E. Belin, F-31400 Toulouse, France 

6 Laboratoire d'Astrophysique de Marseille, Traverse du Siphon, BP 8, 13376 Marseille Cedex 12, France 

7 Physics & Astronomy Department, University of Western Ontario, London, Ontario, Canada, N6A 3K7 



CO 

Abstract. The origin of the strong magnetic fields observed in chemically pe- 
culiar Ap and Bp stars stars has long been debated. The recent discovery of 
magnetic fields in the intermediate mass pre-main sequence Herbig Ae and Be 
stars links them to Ap and Bp stars, providing vital clues about Ap and Bp stars 
and the origin and evolution of magnetic fields in intermediate and high mass 
stars. A detailed study of one young magnetic B star, HD 72106A, is presented. 
This star appears to be in a binary system with an apparently normal Herbig Ae 
star. A maximum longitudinal magnetic field strength of +391 ± 65 G is found in 
HD 72106A, as are strong chemical peculiarities, with photospheric abundances 



of some elements ranging up to lOOx above solar. 



1. Introduction: H Ac-Be and Ap stars 

A few percent of main sequence A and B stars display strong, globally ordered magnetic fields; these 
are the so-called Ap and Bp stars. The origin and evolution of these magnetic fields is currently the 
subject of intensive research. An important avenue of this research is the investigation the progen- 
itors of Ap and Bp stars. Herbig Ae and Be (HAeBe) stars are the pre-main sequence progenitors 
of the main sequence A and B type stars. As such, HAeBe stars have intermediate masses and 
display emission, infrared excess, and tend to be associated with dust and nebulosity. It has been 
long suggested that some HAeBe stars may evolve into magnetic main sequence Ap and Bp stars. 
If this were the case, it was hoped that there would be some distinguishing observable feature in 
HAeBe stars linking them to Ap stars. 

Pioneering observations in 2005 and 2006 successfully detected magnetic fields in HAeBe stars 
for the first time (Wade et al. 2005, Catala et al. 2006). The detected fields appear to display similar 
intensities and geometries when compared to Ap stars (Alecian et al. 2006). The longitudinal field 
strengths detected are on the order of hundreds of gauss, and the geometries are predominantly 
dipolar. Additionally, the incidence of magnetic fields in HAeBe stars is similar to that of the 
magnetic Ap and Bp stars. This strongly supports the idea that magnetic HAeBe stars do in fact 
represent the progenitors of main sequence magnetic Ap and Bp stars. 
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In this paper we investigate one particularly interesting young magnetic B star, HD 72106A, in 
detail. 



2. HD 72106 

HD 72106 is a visual double system, with a 0.8" separation. The primary component has a definitely 
detected magnetic field (Wade et al. 2005), and the secondary is clearly a Herbig Ae/Be (HAeBe) 
star (Vieira et al. 2003). The secondary shows strong, variable emission in Ha, as shown in Fig. [U 

The system has a very large proper motion in a: —5.18 ± 1.08 mas/yr, and in 5: 9.73 ± 1.36 
masy -1 but no relative motion between components is observed (Hartkopf et. al. 1996 & Hipparcos 
data, H0g et al. 2000). The radial velocities of the both components are identical to within ±1 
kms -1 , at 22 kms -1 . As well, the Hipparcos parallax solution for the system places the two stars at 
the same distance: 288(+202/ — 84) pc. This indicates that the two stars are travelling in tandem, 
and strongly suggests that the system is in fact a binary. This is critical because, if the system is 
truly binary, then both components likely have the same age and initial formation conditions. 

Wade et al. (2006) determine, for the primary, an effective temperature of 11000 ± 1000 K, a 
log g of 4.0 ±0.5, and a mass of 2.4 ±0.4 M@. For the secondary they derive an effective temperature 
of 8000 ± 500 K, a log g of 4.25 ± 0.25, and a mass of 1.75 ± 0.25 M®. 

The two stars can be placed on an H-R diagram (Wade et al. 2006), and have consistent ages 
of about 10 Myr, near the zero-age main sequence, when compared with isochrones from Palla & 
Stahler (1993). This is consistent with the components being co-eval. Both the H-R diagram position 
and the association of HD 72106A with the HAeBe secondary provide evidence that the primary is 
a pre-main sequence or very early zero-age main sequence star. The H-R diagram of both stars is 
shown in Fig. [2 
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Figure 1: Observed Ha profiles for HD 72106B (the apparently non-magnetic secondary) from 11 
Jan. 2006 (dark/black) and 11 Feb. 2006 (light /green) . Strong and variable emission is observed, 
indicating the presence of an emitting circumstellar envelope. 
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Figure 2: The pre-main sequence H-R diagram of the HD 72106 system. The solid lines are the model 
evolutionary tracks of Palla and Stahler (1993), with masses as labeled, and the dot-dashed lines 
are isochrones, with ages labeled. The birthline (which assumes an accretion rate of 10 _5 MQyr _1 ) 
and zero age main sequence (ZAMS) are shown as dashed lines. 



3. Spectrum Reconstruction 

Ten Stokes V spectra of HD 72106 were obtained between February 2005 and February 2006 with 
ESPaDOnS, the high-resolution echelle spectropolarimeter mounted at CFHT. The spectra cover a 
wavelength range from 3700 Ato just past 1 /xm. Problematically, the ESPaDOnS pinhole is 1.6" 
in diameter, whereas the two components of HD 72106 are separated by 0.8". With good seeing 
and careful guiding it is possible to isolate individual components. This was accomplished on one 
occasion. However, the majority of our spectra are of combined light from the system. Observations 
were conducted using atmospheric dispersion correction, with careful attention to keeping the system 
centroid centred in the pinhole. 

In order to isolate the spectrum of the primary star, we subtracted the individual spectrum of 
the secondary (obtained on a night of excellent seeing) from the combined spectrum, weighted by 
its luminosity contribution. Assuming the stars are located at the same distance, their apparent 
magnitudes, suitably bolometric-corrected, yield their relative luminosities. The relative apparent 
magnitudes of the components of HD 72106 are well known (eg. Fabricius & Makarov 2000), and the 
Hipparcos parallax solution puts both components at the same distance. The subtraction process 
also assumes that the spectrum of the secondary is non-variable, which is a reasonable assumption 
for metallic lines of a non-magnetic late A star (excluding Ha). Thus, with the relative apparent 
magnitude, we can perform the subtraction effectively. 

To check this method, Balmer lines (excluding Ha), as well as several deep, apparently non- 
variable metallic lines, were compared in the recovered spectra and in the individual observed 
spectrum of the primary. A very good agreement was found. Presented in Fig. [3] is a typical ob- 
served spectrum of the combined system, compared with a spectrum of just the secondary, and the 
reconstructed spectrum of the primary, resulting from the subtraction process, compared with a 
spectrum of just the primary. 
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Figure 3: The top frame shows a sample of the spectrum from the secondary (in black) and an 
observed spectrum of the combined light from the system (in gray). The bottom frame compares the 
spectrum reconstructed spectrum of the primary (in gray), resulting from the weighted subtraction 
process, with an observation of just the primary (in black). A very good correspondence is obtained 
(although some differences are expected due to the line profile variability of the primary). 



4. Magnetic Field Measurements 

Least Squares Deconvolution (LSD) (Donati et al. 1997) was performed on the reconstructed spec- 
tra of the primary. LSD is a method of effectively averaging over many lines in a spectrum in order 
to improve the signal to noise ratio of the line profile. A line mask for a 10000 K Bp star, and 
spectra with the contribution of the secondary removed, were used. Significant Stokes V signatures 
are detected at all observed phases. Longitudinal magnetic field measurements were obtained from 
the LSD Stokes V profile for each observation. The largest longitudinal field found was +391 ± 65 
G. Stokes V profiles with larger amplitudes were observed, but with crossover signatures, corre- 
sponding to weaker longitudinal fields. The Stokes I and V LSD profiles corresponding to the 391 
G measurement are shown in Fig. 01 and are labeled 0.4015 in the phase column. 



5. Period Searching 

Clear variability, presumably due to rotation, can be seen in the spectrum of HD 72106A, with an 
apparent period around two days. To investigate this further we first searched our measured longi- 
tudinal magnetic field data for a period. This was done using a modified Lomb-Scargle technique, 
fitting sinusoids through the data to produce a periodogram, and looking for a minimum in x 2 - 
Unfortunately, due to the small number of longitudinal field observations, a unique solution could 
not be found. 

In an attempt to improve our period determination, a more sophisticated technique was used. 
This method is still based on fitting a sinusoid through a time series of data points to produce a 
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periodogram. However, now the set of measurements for one pixel in the LSD profile is used as the 
time series of data points. The periodogram for one data point may be heavily affected by noise, 
but if one averages the periodograms of all the pixels in the LSD line profile, producing an average 
periodogram, the effects of noise are substantially reduced and the true period found. To verify the 
technique, periods were determined for a number of known Ap stars (Auriere et al., in preparation), 
as well as synthetic test spectra, with good results. 

While this improved the situation significantly by reducing the number of possible periods, it 
still did not produce a unique period solution for HD 72106A. Examining phased LSD profiles by 
eye allowed us to reject several more periods as un-physical, but a few possibilities still remain 
around 1.8 and 1.9 days. A set of LSD profiles phased with one of the better candidate periods is 
presented in Fig. HI 

6. Spectrum Synthesis 

Preliminary spectrum modeling has been performed for both the primary and the secondary, for a 
single phase only, using observations of the individual components. The Zeeman2 spectrum synthesis 
code was used, with ATLAS9 model atmospheres corresponding to the temperatures and gravities 
deduced by Wade et al. (2006). In the primary, a dipole magnetic field of 1 kG was assumed, 
oriented along the line of sight, with no microturbulence. A homogeneous abundance distribution 
both vertically and horizontally was assumed. In the secondary, a homogeneous vertical and hori- 
zontal abundance distribution was used, with no magnetic field. Microturbulence was left as a free 
parameter in the fitting procedure. 

We find that the primary displays clear chemical peculiarities. The abundance pattern is typical 
of those of Ap/Bp stars, although the abundance enhancements are remarkably strong. In particular, 
strong overabundances of Cr, Fe, Si, and Ti are found. In contrast, the secondary displays solar 
abundance for all elements detected. 

Clear structure and variability can be seen in the line profile of the primary, due to abundance 
spots on the surface of the star. These appear most clearly in lines of Fe and Cr, but can also be 
observed in some Ti and Si lines. Complex, variable line profiles of this type are typical of Ap stars, 
and are explained by the presence of abundance non- uniformities in the stellar photosphere. The 
secondary displays smooth line profiles with no visible asymmetries (although the S/N of this star 
is relatively poor, and additional observations are required to confirm this). 

A sample spectrum of the primary, with a best fit synthetic spectrum, is shown in Fig. [5j Best fit 
abundances and v sini for the primary are shown in Table [TJ For the secondary, best fit abundances, 
vsinz, and microturbulence can also be seen in Table [U with sample observed and synthetic spectra 
in Fig. [6] . 

7. Conclusions 

These results demonstrate that strong chemical peculiarities, of the type seen in Ap/Bp stars, 
together with magnetic fields, can exist in late B-type stars at the pre-main sequence or very early 
zero-age main sequence phase. A maximum longitudinal magnetic field of +391 ± 65 G was found 
in HD 72106A. Overabundances of up to 2 dex above solar were also found, along with a complex 
absorption line structure, suggesting an inhomogeneous surface abundance distribution of elements. 
These properties are all characteristic of Ap stars. In contrast, HD 72106B which presumably formed 
at the same point in time, under the same conditions as the primary, displays solar abundances and 
no magnetic field. Thus the HD 72106 system provides a unique link between Ap stars and their 
pre-main sequence progenitors. 

Further circular polarisation observations of the HD 72106 system are currently planned. This 
will allow us to better sample the spectroscopic and polarimetric cycle of the primary to determine 



6 



v/i 



phasG 




0.0189 



0.2178 



0.2383 



0.2551 



0.2942 



0.4015 



0.4206 



0.721 1 



9202 



-50 50 100 -50 50 100 
Velocity (km/s) Velocity (km/s) 



Figure 4: The Stokes I and V profiles obtained for HD 72106A, phased according to a 1.953 day 
period. The Stokes V profiles, detected significantly at all phases, have been multiplied by 20 times 
relative to the I profiles, and shifted vertically for clarity. The profiles are labeled by phase, with an 
arbitrary zero point. Clear variations of the Stokes I profile are apparent. 



a unique rotational period. With additional observations we can model the magnetic field geometry 
and map the surface distributions of elements. Additional observations will also allow us to test the 
stability of the spectrum of the secondary and further test the spectrum reconstruction technique. 
With these new observations we will be able to more completely describe this unique and remarkable 
system. 
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Element 


Primary 


Secondary 


Al 


-0.2 ±0.4 




Si 


+1.0 ± 0.2 




Ca 




-0.4 + 0.4 


Ti 


+1.0 ±0.3 


-0.1 + 0.4 


Cr 


+1.9 + 0.1 


+0.1 + 0.3 


Fc 


+0.9 + 0.1 


-0.1 + 0.3 


Ba 




-0.3 + 0.5 


v sini 

e 


41 + 3 kms" 1 
~ kms -1 


50 + 4 kms" 1 
2.5 + 0.5 kms" 1 



Table 1: Preliminary best-fit abundances, wsini, and microturbulence values for HD 72106A and 
HD 72106B. Abundances are based on the 4500-4700A region. 
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Figure 5: Observed (black) and best fit synthetic (gray) spectra for HD 72106A. Major contributers 
to each feature have been labeled. 
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Figure 6: Observed (black) and best fit synthetic (gray) spectra for HD 72106B. Major contributers 
to each feature have been labeled. 
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